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ABSTRACT 

Aims. We try to identify the origins of field O-stars in the nearest 2 to 3 kpc around the Sun using the best presently available kinematic data 
on O-stars and on young open clusters. We investigate the question if the present-day data are consistent with the assumption that O-stars have 
formed in groups (clusters, associations), or in isolation. 

Methods. We apply the epicycle theory for back-tracing the orbits of O-type stars and of candidate parent open clusters. 
Results. From the 370 O-stars in the "Galactic O star catalog v 2.0" (GOSV2) we have investigated 93 stars classified as field, and found the 
origin for 73 of them in 48 open clusters younger than 30 Myrs. Only for 32 stars or about 9% of all O-stars from this catalogue, the question 
of their origin in groups is not solved; some of them may have originated in isolation or may have disintegrated the group in which they 
formed. Fifty percent of the young open clusters (age < 30 Myr) in the "Catalogue of Open Cluster Data" (COCD) have O-stars as members, 
or have ejected at least one O-star in the first 10 Myrs of their life, or both. During this period the average mass loss from open clusters by 
ejecting O-stars is found to be 3 to 5 M Q per Myr. We prove that f Pup had its origin in the open cluster Trumpler 10 which it left about 2.5 
Myrs ago, and that its present-day distance is 300 pc (compared to 440 pc before). The revised distance implies a significant revision of the 
stellar parameters (a radius of 14 R Q , a mass of 22.5 M Q , and a luminosity of log L/L G of 5.74) i.e, ( Pup is closer, less massive, and less lumi- 
nous than previously thought. Our findings provide independent estimates of the present-day distances and absolute magnitudes of field O-stars. 
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1. Introduction 



Do all O-stars form in groups (clusters, associations) as 
is commonly believed or is the formation of O-stars in iso- 
lation possible? This long-standing question can only be an- 
swered, when the birth-places of all O-stars will be discov- 
ered. A revi ew of the situati on is given in the in t roduct i on by 
iGiesI d 1987b and recently in IZinnecker & Yorkd d2007l) . iGiesI 
dl987l) compilied a catalogue of 195 O-stars which he used to 
infer the first solid statistics about runaway and field O-stars. 
R ecently, a new catalogue of G alactic O-stars (GOS V1 version 



1, M afz-Apellaniz et all d2004l) : GOSV2 version 2, ISota et al. 
d2007ln was published. Comprising 370 O-stars, the catalogue 



allows to re-address the statistics of O-star birth-places. In par- 
ticular, the GOSV2 catalogue contains a subset of 105 O-stars 
called^eW, which simply means that they cannot be identified 
as present or former members of recognised groups. Such "iso- 
lated O-stars" are of essential interest to decide the question if 
"isolated" massive star formation is possible or not. 

Because of the relatively short lifetime (a few million years) 
near the main sequence, the orbit of an O-type star in the wider 
solar neighbourhood can, in principle, be followed all the way 



back in time to the onset of its hydrogen-burning stage. This 
means that the location of its parent star forming cloud can be 
determined. Investigating the area around these parent clouds, 
one may find other young objects there, e.g. young star clusters 
or OB-associations. 

During the last decade, after the results of the Hipparcos 
mission became available, nearby OB -associations have bee n 
investigated in considerable detail dde Zeeuw et al. Il999t) . 



However, not all OB-stars have been found living in associa- 
tions, some are far way from presently-known stellar groups 
on the sky. 

Using the data from ES A's Hipparcos mission, 



Hoo gerwerf et al.l (2001) back-traced the orbits of 56 OB-type 
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runaway stars and nine compact objects with distances less 
than 700 pc. They found that at least 21 objects of their sample 
could be linked back to nearby associations and young open 
clusters. The authors state that the remaining objects may have 
originated from distances farther away than 700 pc, where the 
knowledge of parent groups is poor. 

Another line of argumen t has been followed in two papers 
by Ide Wit et alJ (12004, 2005). In their first paper they investi- 
gate the origin of 4 3 O-type field stars from the O-star cata- 
logue bv lGiesI dl9871) . The authors search the area around these 
stars for stellar groups in the near-infrared which are possibly 
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hidden in the optical. In their second paper (Ide Wit et al . 2005) 
they investigate the same sample kinematically. They excluded 
as field stars the runaway stars, i.e. those with spatial y eloc- 
ities above the limit of 40 km s _1 set by iBlaauwl d 19611) and 
those at distances larger than 250 pc from the Galactic plane. 
Combining the results of both papers they claim that not more 
than 4 + 2% of all O-stars in Gies' catalogue can be called 
genuine field stars. 

The argument by Ide Wit et al.l (2004, 2005) is a rather in- 
direct one, they are excluding stars as field without being able 
to retrace their origin. This is exactly the point where we start 
our present study. Only if one succeeds to retrace an O-star to a 
parent group within its past lifetime one can say with certainty 
that this O-star has originated in a group. Proving or disproving 
this point is not an easy task given our incomplete knowledge 
of possible birth-places in the wider neighbourhood of the Sun 
and the uncertainties of the six-dimensional phase space coor- 
dinates (position and motion) of candidate stars and candidate 
clusters and/or associations. 

In this paper we are testing the hypothesis that O-stars, the 
origin of which is hitherto unknown, may have been ejected 
from young open clusters (or their protoclusters) during or after 
the star formation period in the (parental) cluster. For this pur- 
pose, we follow the path of stars and clusters back in time in the 
Galactic potential. In the next section we present the underlying 
observations, then we describe the method and its application. 
Section |4] is a presentation and a discussion of the results. In 
Section [5] we consider a few selected examples of stars with 
the adopted solutions, whereas in Section [6] we briefly discuss 
the stars for which we did not find an acceptable solution. A 
summary concludes the paper. 



2. Observational material 



For the back-tracing of stellar and cluster orbits we are us- 
ing the most homogeneous and accurate data of all 6 param- 
eters of space phase available at present. We took the posi- 
tions and proper motions fro m the recently completed PPMX 
catalogue (Ros er et al.lE o08). and the radial velocities from 
the CRVAD-2 (Kh archenko et al.ll2007l) . The major sources of 
specific informati on on open clusters and O -type stars were 
the ca talogues by iKharchenko et al. I J2005allbh and ISota et al 
(12007b . respectively. 



2.1. Open clusters 

The C atalogue of Open Clust er Data (COCD) and its 
Extension (K harchenko et al. 2005a b) includes 641 open clus- 
ters and 9 cluster-like associ ations identified in the ASCC-2.5 
catalogue dKharchenkol200ll) . For each cluster the membership 
was d etermined using spatial, kinematic, and photometric cri- 
teria ( IKharchenko et al.l 120041) . A homogeneous set of cluster 
parameters was derived by applying a uniform technique. The 
nine associations are included in the COCD because of their 
compact appearance on the sky and small dispersion in proper 
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Fig. 1. The distribution of the rms errors of proper motions 
H a cos 6 (left), (middle), and of radial velocities RV (right). 
Upper panels show the distributions for the 161 young open 
clusters, lower panels for the 93 field O-stars. The peak at 10 
km/sec in ctrv of O-stars represents stars for which no infor- 
mation on the rms error of the radial velocity is given in the 
literature. 



motion space, so that they can be treated via our standard mem- 
bership selection procedure. The COCD contains the celestial 
position of a cluster, its distance to the Sun, reddening, age, 
angular size, proper motions, and, if available, radial velocity. 
Recently, th e parameter set was supplemented by tidal radii 
and masses dPiskunov et al. 2008a). The completeness of the 
cluster sample is mainly defined by the limiting magnitude of 
V ~ 11.5 of the ASCC-2.5. Therefore, even nearby embedded 
clusters can been missing in the sample if their members were 
fainter than V ~ 1 1 .5 in the optical. Nevertheless, for "clas- 
sical" open clusters, i.e. when the bulk of the placental matter 
is removed and clusters become visible in the optical spectral 
range, the sa mple was found to be c omplete up to a distance of 
about 850 pc dPiskunov et al.ll2006l) . 

In the context of this paper, however, we are not interested 
in the full sample, but in a sub-sample of young clusters. As 
potential candidates for parent groups we considered 161 clus- 
ters (including 9 associations) younger than 30 Myr, and hav- 
ing radial velocities measured. Since absolutely bright stars are 
still present i n these clusters, this subset is volume limited to 
about 2 kpc ( Piskunov et al. 2006b. o r to a distance modulus 
(V - M v ) « 13 dSchilbach et al.ll2006b when extinction is taken 
into account. Except for proper motions, the cluster data for 
each cluster (i.e., coordinates of the cluster centre, distance, and 
radial velocity) were taken from the COCD. The mean proper 
motions were recomputed fr om the PPMX data and the mem- 
bership information given in Kharchenko et al.l ( 2004 ). In the 
upper panels of Fig. Q] we show histograms of the mean errors 
of the kinematic data for this sub-sample of young clusters. 



2.2. O-type stars 
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The sample of O-type st ars for this pape r was taken from 



the second version GOSV2 (Sota et al 



20071) o 



of the "Galactic 



O star catalog" by Ma fz-Apellaniz et a l. (2004), with 370 en- 
tries. The catalogue is expected to be complete for O-stars 
brighter than V — 8 but it includes many fainter stars, too. For 
each star, the catalogue delivers spectral classification, photo- 
metric and astrometric data, and further information such as 
multiplicity and membership in known associations. Among 
370 stars, 105 stars are classified as field stars or field runaways 
from unknown parent groups. These are the stars we were in- 
terested in. 

Again, we took the positions and proper motions from the 
PPMX catalogue. Doing so, we not only benefitted from ac- 
curate proper motions but also kept proper motions of clus- 
ters and O-stars on the same system. Radial velocities came 
from CRVAD-2. Though CRVAD-2 presents data for 55000 
stars, radial velocities are available for only 93 O-stars of our 
list. Unfortunately, for 41 stars, no information on the nns er- 
rors of radial velocities is available, so we assumed them to be 
+ 10 km/s. The histograms of the mean errors of kinematic data 
for the final sample of O-stars are shown in the lower panels of 
Fig. 03 As expected, the accuracy of the data for clusters is, on 
average, higher than for O-stars. For clusters, the median of the 
nns errors of fi a cos 5, fig and RV are 0.5 mas/yr, 0.4 mas/yr and 
3.4 km/s, respectively. For O-stars, the corresponding numbers 
are 1.2 mas/yr, 1.0 mas/yr and 7.5 km/s. 

The distances of O-sta rs are m ore of a problem. The 
GOSV2 gives Hipparcos dESAl 19971) p aral laxes for all s t ars for 
which these are available. Very recently. Ivan Leeuwenl ([2007) 
published the new reduction of the Hipparcos observations. 
Due to a sophistica ted modelling of the satellite's attitude, 
van Leeuwenl (120071) could considerably improve the trigono- 



metric parallaxes of stars brighter than about 8th visual magni- 
tude. However, only five O-stars of our sample had Hipparcos 
parallaxes with an accuracy better than 30%. 

In order to get distances for all stars, there was no other 
chance than using the methods of distance estimates based 
on the spectroscopic and photometric data. To derive spectro- 
scopic distances d sp from the well known relation 



logd sp = 0.2[V -M v + 5- 3.1((B -V) — (B- V) )], 



(1) 



we took the observed B and V magnitudes and the spectral 
classification from the GOSV2 catalogue, and converted spec- 
tral type into ab solute magnitude My and (B-V)o according to 
Schmidt-Kalerl dl982r) . However, the spectral classification of 
O-stars is neither straightforward nor unambiguous. For ex- 
ample, for a relatively bright star, HD 135240 (V = 5.08), 
one finds a spectral type and luminosity class of O 7.5 III in 
the GOSV2 catalogue and O 8.5 V in SIMBAD. This differ- 
ence in spectral classification leads to an uncertainty of about 
1 mag in distance modulus. The second source of uncertainty 
arises from the My-spectral-type calibration which can intro - 
duce systematic effects of up to 1 mag (e.g., see lWalbornl2002l) . 
Moreover, the calibration of O-stars shows a l arge scatter which 
may be intrinsic to the star s themselves (e.g., Conti et all 1983 



Garmany & Stencellll992l) . We conclude that the spectroscopic 
distance moduli may be uncertain by up to 2 mag, and the dis- 



tance of a field O-star is the most inaccurate input parameter in 
the back-tracing of stellar and cluster orbits. 

3. The back-tracing method 



3.1. The epicycle approach 



For the re-tracing of the stars and open cl usters we fol- 
lowed the approach used by iFuchs et al. (2006) who a dopted 
the epicyclic equations of motion as given by Eindbladl (fl959): 



m = 
n(t) = 

((f) = 

u(t) = 

v(t) = 
w(t) 



i/(0) k(0) 
2B sin(/cf) 

7 (0)+2A[£(0)+^]*- 

ZD 

-—i/(0) sin(/tf) + — h(0)[1 - cos(tf)], 
Bk k2 



(2) 



w(0) 

((0) cos(vf) + — — - sin(vf), 



m(0) cos(/tf) + — v(Q) sin(Kf), 
2B 

2B 

m(0) sin(*t) + i/(0) cos(/cf), 

K 

w(0) cos(vf) - £(0)vsin(vr) 



(3) 



where A and B are Oort's constants, Qo is the angular ve- 
locity of the Galactic rotation of the local standard of rest 
(LSR), k and v are the epicycle frequency and the vertical os- 
cillation frequency, respectively. We assumed a flat rotation 
curve A - —B = Qo/2, where Q.q = 25.9 km/s/kpc at the 
Galactocentric radi us of the LSR ro _ 8 .5 kpc. We also adopted 
the same values as IFuchs et al. ( 2006b for the parameters k — 
39.0 km/s/kpc and v — 74 km/s/kpc. 

Eqs. (2) and (3) describe the motion of a particle in a 
non-inertial coordinate system centred at a fiducial point at a 
Galactocentric radius ro (at Z=0) from the Galactic centre, for 
which the transformation from the Cartesian Galactic coordi- 
nates X, Y, Z into j], £ is given by 



( = r -r, 

j] = ro x arctan(T/r), 

( = z, 



(4) 



with r = -y/(ro - X) 2 + Y 2 . Similarly, the velocity components 
U, V, W of the peculiar space velocity (after correcting for solar 
motion and Galactic rotation) are transformed into u, v, w via 



r -X Y 
u = U— V-, 

r r 

Y ro-X 

v = U- + V- , 

r r 

w = W. 



(5) 



As we are only interested in the relative location and ve- 
locity of a star with respect to a candidate parent cluster, 
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we chose r close to the Galactocentric radius of the can- 
didate cluster in each case, and the flat rotation curve gave 
D.0 = 220/ro km/s/kpc. The epicycle and the oscillation fre- 
quencies have been assumed to be constant. This extends the 
validity of the approach, as the initial rj are small for star 
and cluster, even at larger distances from the Sun, provided that 
r s tar - ro « ro, and v « ro £2o- We checked that the require- 
ments were fulfilled for our solutions. 

3.2. Variation of the initial conditions and selection of 
the solution 



As a starting point of the backward computations, we took 
the positions, proper motions, radial velocities of stars and 
clusters at their nominal values from the sources described in 
Sec. [2] However, taking into account the quality of the input 
data, we allowed variations of initial conditions for O-stars in 
eqs. (2) and (3) within the given error budget. We varied the 
nominal values of their proper motions and radial velocities 
between -2.5 <x and +2.5 cr in 0.5 <x steps, and the distance 
moduli between [(V - My) - 2] and [(V - My) + 2] in steps of 
0.02 mag. Initial conditions for the clusters were not varied. 

The orbits were traced back in time over 1 1 Myr with a 
step of 0.05 million years. At each time step, the relative dis- 
tance between each star and each cluster was determined. For 
the rest of the paper we adopt the following terminology. We 
speak of an encounter, if, going backward in time, the distance 
between star and cluster (centre) was less than 10 pc. This is 
called an acceptable solution. The encounter time, t enc , is the 
time before present when an encounter occurred, and the en- 
counter distance, d enc , is the distance between star and cluster 
centre at t enc . Note, however, that the actual physical process is 
an ejection of a star from its cluster counterpart. This ejection 
occurred at time t enc , when the star was at a distance d enc from 
the cluster centre. 

For each acceptable solution, we computed a probability 
Pkin describing how well the proper motions and radial veloc- 
ities (actually used in the solution) suit the nominal values of 
the kinematic parameters of given star and cluster: 



1 

Pkin = exp<-- 



Afi a , s \ I Ap StS 



ARV, 



(6) 



Here Afi a<s , Afj,g s , and ARV S are differences between the applied 
and the nominal values of the corresponding velocity compo- 
nents of the star, and (E vc ,s,c) 2 = o"? c s + o- 2 vc c where cr vc ;S and 
cr vc , c are the rms errors of the corresponding velocity compo- 
nent of the star and cluster under consideration. If A is smaller 
than E in each of three velocity components, then the corre- 
sponding solutions has a kinematic probability pki„ > 0.606. 

This procedure yielded a three-parameter set i.e., encounter 
time t enc , present-day (or dynamic) distance ddyn of an O-star 
from the Sun, and kinematic probability pki„ of acceptable so- 
lutions for the encounter distance d mc . In the next step, we se- 
lected 10% of solutions with the largest p&n and gave the high- 
est priority to the solution providing the smallest d enc . If, for a 
given O-star there were solutions hitting more than one clus- 
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Fig. 2. Dynamic distance versus back-tracing time for 
HD 75222 (left panel) and HD 201345 (right panel). Dots 
indicate all acceptable solutions, stars show the best solu- 
tions for an individual star-cluster combination. Vertical lines 
mark different distance estimates for these stars: the solid 



fine is for a spectro-photometric distance from iNeckel et al 
(1980), the dotted and dashed lines ar e for spectroscopi c dis- 
tances based on spectral types fr om ISota et alj (12007) and 
My-spectral-type ca l ibratio ns from ISchmidt-Kalen (fl 982) and 
Garmanv & Stencell dl992l) . respectively. The dashed-dotted 
line is a dist a nce e stimate taking the spect ral type from 
Houk & Swift ( 19991) and calibrations from ISchmidt-Kalei 
dl982h . 



ter, we took into account the available information on distance 
estimates for this star (e.g., spectroscopic distance, spectro- 
photometric parallax, Hipparcos parallax) for making the de- 
cision. 

In Fig. [2] we give two examples of such sets of solutions. 
The left panel of Fig. |2]shows the case of HD 75222, the right 
panel is for HD 201345. For HD 75222 we obtain three sets 
of acceptable solutions depending on the present-day distance 
ddyn of the star from the Sun. Nevertheless, the solution includ- 
ing vdBergh-Hagen 56 as a cluster which the star encountered 
about 3 Myr ago can be rejected due to a bad compatibility 
with spectroscopic and spectro-photometric distance estimates. 
The solution with Trumpler 15 fits a distance estima t e base d on 
My-spectral-type calibrations from ISchmidt-Kalei (1982) but 
needs strong variations in all three kinematic parameters, so 
even for the best solution the kinematic probability is smaller 
than 0.2. On the other hand, the solution with Collinder 205 
coincides better with the other distance estimates available for 
HD 75222, and it is quite a stable one. We conclude: assum- 
ing a present-day distance dd yn - 1735 pc, HD 75222 was 
ejected from the young cluster Collinder 205 at d enc * 1 pc 
about 6.6 Myr ago. The kinematic probability of this solution 
is 0.96. 

The case of HD 201345 is less clear. There are four sets 
of acceptable solutions. Again, encounters with NGC 6913 and 
Berkeley 86 can be rejected since they assume a present-day 
distance of the star of about 1.2 kpc which is too small with 
respect to the spectroscopic distance estimates. The solutions 
with NGC 687 1 and Cyg OB2 as counterparts assume both a 
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(B-V)o, mag 

Fig. 3. The colour-magnitude diagram ((2? - V)o, My) of the ac- 
ceptable solutions for HD 75222 and HD 201345. The dashed 
lines show (My — 3.1(2? — V)o) lines based on eq. (1) and dy- 
namical distances dd yn derived from the back-tracing proce- 
dure for HD 75222 and HD 201345. The dotted line stands 
for the ZAMS, whereas solid lines indicate isochrones of 
Collinder 205 and Cyg OB2 found to be the hosts of HD 75222 
and HD 201345, respectively. Small circles mark the abso- 
lute magnitude adopted for HD 75222 (My = -5.6) and 
HD 201345 (My = -4.0). 



present-day distance of dd vn = 1740 pc for HD 201345 and 
fit much better the spectroscopic distance estimates. With a 
kinematic probability /%„ = 0.99 the star was ejected from 
Cyg OB2 at d enc = 1.5 pc about 5.8 Myr ago. For the solu- 
tion with NGC 6871 we obtained d enc = 1 pc, t enc = 6.5 Myr, 
Pkin = 0.85. Since the kinematic probability of the solution with 
Cyg OB2 is higher, we selected Cyg OB2 as the most proba- 
ble host of HD 201345. However, the solution with NGC 6871 
cannot be rejected completely. We discuss this case below in 
Sec. [5] which is devoted to the results on individual stars. 

For each star with a present-day distance dd w determined, 
we can estimate its absolute magnitude My using eq. (1) which 
directly provides a relation (My - 3.1(2? - V)o). Assuming that 
the adopted extinction law is correct and the intrinsic colour 
(2? - V)o of an O-star is between -0.33 mag and -0.13 mag, 
we obtain a maximum variation of about 0.6 mag for its ab- 
solute magnitude My. The limit of -0.33 mag is defined by 
the location of the ZAMS, whereas the limit - 0.13 mag corre- 
spond s to spectral types from B7 V to B3 I (see lSchmidt-Kalei 
and it was chosen to be on safe ground not to exclude 
O-stars. Assuming further that the star is "genetically" related 
to its counterpart-cluster, we chose the crossing point between 
the line (My - 3.1(2? - V)o) for the star and the isochrone 
corresponding to the cluster age to be the absolute magnitude 
My of the star. For illustration, we show the corresponding 
((2? - V)o,M v ) plot in Fig. 3 for HD 75222 and HD 201345. 
The adopted absolute magnitudes My are -5.6 and -4.0 for 
HD 75222 and HD 201345, respectively. 
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Fig. 4. Distribution of off-sets introduced in proper motions 
yU a cos<5 (left), (middle), and in radial velocities RV (right) 
of the O-stars getting acceptable solutions. 

4. Results 



For 73 out of 93 O-stars considered, we found acceptable 
solutions indicating that the present-day data are consistent 
with the assumption that these O-stars had encountered (ac- 
tually, are ejected from) young open clusters during the past 10 
Myr. The essential results of the back-tracing calculations are 
compiled in Table [2] For each of the 73 O-stars we give: HD 
identification (column 1), its spectral classi fication taken from 



the GO SV2 catalogue (2), name and age (IKharchenko et al 



2005alTbh of the probable counterpart cluster (3,4), the dynam- 
ical distance of the star dd yn (5). Column 6 gives the time t enc , 
i.e. the time before present when the star was ejected. Columns 
7 and 8 are the distance d enc and the relative velocity AVel enc 
between the star and the cluster at t enc . Column 9 contains the 
star's absolute magnitude and its upper and lower limits esti- 
mated from dd w (see the end of Sec. l3.2t . and column 10 gives 
the kinematic probability of the solution. 

Except in one case, the kinematic probability p^ m is always 
larger than 0.5. This means that in units of rms errors, only 
small variations of the kinematic parameters were needed to 
get an acceptable solution for the majority of stars. The dis- 
tributions of off-sets in the sense (used parameter - nominal 
parameter) are shown in Fig. [4] for proper motions and radial 
velocities of the 73 stars with acceptable solutions. The dis- 
tributions do not indicate any anomalies, so we conclude that 
the results of the back-tracing calculations can be used to un- 
derstand where the field O-stars came from and explain their 
present location. 

Of 161 young open clusters in our sample, only 48 clus- 
ters were "hosts" of field O-stars. Only two clusters are older 
than 20 Myr. Twelve clusters ejected two O-stars each, two 
clusters (Cyg OB2, and ASCC 8) had three ejections in the 
past, and three clusters (Trumpler 14, Loden 821 and ASCC 9) 
even four ejections. Three O-stars turned out to be members 
of the newly detected op en clusters ASCC 45 and ASCC 79 
( Kharchenko et al . 2005b). Since their origin seems to be clear, 
they are excluded from the statistics presented below. 

In Fig. [5] we show the ages of the clusters versus encounter 
time t enc . Close to the bisector in this figure, one expects to find 
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Fig. 5. Cluster age f oc ; versus encounter (ejection) time t enc . 
Dots stand for the 70 O-stars with acceptable solutions from 
the back-tracing procedure. Crosses mark 23 O-stars which 
were possible cluster members at the beginning of their life 
but, later, were ejected from the parent clusters. The solid line 
tenc - 0.5 x t oc i separates them from the other O-stars whereas 
the dotted line is the bisector. The red open circles indicate O- 
stars with velocities larger than 80 km/s with respect to their 
cluster counterparts at the moment of the closest approach 
(ejection). 



O-stars which left the cluster birth-places just before, during 
or just after the time the clusters were forming. In this sense 
they were members of protoclusters but did not become ac- 
tual cluster members; they were formed in the same regions 
where the corresponding clusters originated. On the other hand, 
there is a number of O-stars with t enc significantly smaller than 
the ages t oc \ of their counterparts. For these stars, we assume 
that they were ejected from the already formed cluster at an 
early stage of the cluster's life, either due to internal evolu- 
tionary processes in the cluster itself, due to bi nary evolution 



or due to external dist urbing forces (see e.g. iRroupa & Boily 



2002; Zinneckerl 12003 ). We have marked these candidates by 
crosses in Fig. when t enc < 0.5 x t oc \. This criterion is re- 
lated to the accuracy of the age estimates of open clusters in 
our sample which is fou nd to be about cr\ ogt - 0.20. ..0.25 
(Kha rchenko et al.ll2005al) . We remark that this determination 



of clu ster ages is based on the Padova isochrones (IGirardi et al 
20021) which have a lower limit at log t oc i = 6.6. Therefore, the 



ages of the youngest clusters can be somewhat overestimated. 
Taking this into account, the portion of ejected O-stars from al- 
ready formed clusters should be about 25 - 35%. The majority 
of O-stars, however, was ejected from the star formation region 
during the protocluster phase. The fact, that we do presently 
observe the outcome of this latter scenario as an "isolated" O- 
star and a surviving open cluster, may indicate that cluster dis- 
ruption by O-st ars in early stages is p ossibly less effective than 
assumed by e.g. Lada & Ladal (120031) . 

From the study of the initial m ass function of Galactic open 
clusters, Pisk unov et al. (2008b) found that a typical cluster 
loses about 60 - 80% of its initial mass during the first 260 
Myrs of its evolution. T he average mass loss rate determined by 
Piskunov et al.l d2008bl) ranges from 3 to 14 M /Myr which in- 
cludes mass loss due to stellar and dynamical evolution. From 




n . , . n 



AVel, km/s 

Fig. 6. Distribution of relative velocities of O-stars with respect 
to their cluster counterparts at the moment of the closest ap- 
proach (ejection). The filled histogram shows former cluster 
members ejected from their hosts. 



the number of the ejected O-stars, the number of the parent 
clusters, and the distribution of their encounter times t enc , we 
can roughly estimate an average mass loss of a cluster caused 
by the ejection of O-stars alone. Assuming a typical mass of an 
O-star of 20 M , we obtain the average mass loss rate of a clus- 
ter due to ejected O-stars from about 5 M Q /Myr if t enc < 5 Myr 
to 3 M /Myr if 5 Myr< t enc < 1 1 Myr. Though the statistics 
is rather poor, the resu lt agrees well with the estimation by 
Pisk unov et al. (2008b) and underlines the importance of the 
contribution of ejected O-stars to the general mass loss of open 
clusters in the first 10 Myr of their life. 

At the moment of ejection, the relative velocities of for- 
mer cluster members with respect to their parent clusters are 
rather moderate, with median at AVel « 45 km/s. This does not 
change significantly if we consider all 70 O-stars. In Fig [6] we 
show the histogram of relative velocities AVel for all O-stars 
with acceptable solutions and for the 23 stars probably hav- 
ing been cluster members in the past. Of five stars with AVel 
larger than 80 km/s, there is only one which was a former clus- 
ter member. The other four belonged probably to protoclusters 
and were ejected from the regions during cluster formation. 

As described in Sec. 4, our back-tracing procedure allows 
to determine the present-day, or dynamical distances dd yn of 
O-stars. This gives us estimates of their absolute magnitudes 
Mv,dy n - In Fig.[7jwe show the distribution of diff erences AMy = 
My st — Myd vn to t he stan dard calibrations from Schmidt-Kalerl 



(11982b . IWalbornl dl972l) . and iGarmanv & Stencell ( 1 1992b . At 
first glance this picture is intriguing, but remember that the cal- 
ibrations of absolute magnitudes of O-stars date back to the 
pre-Hipparcos era. After Hipparcos it has become possible to 
start a re-calibration, admittedly more for early B-type stars 
than for the O-sta rs themse l ves. Utilizing the measurements 
from Hipparcos, iKaltcheval (11999b found My, s t ~ M v>Hlpp = 
-0.85 +0.12 for 44 B0-B3 III, IV stars. This number coincides 
well with our findings for the median of AMy = -0.84, -0.72 
or -0.36 depending on the different calibrations (see Fig. 17). As 
we will show below, ( Pup, the only O-star with a highly signif- 
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Fig. 7. Distribution of differences in the absolute magni- 
tudes of field O-stars derived from our dynamical dd yn dis- 
tances versus the absolute m a gnitud e s according to th e cal- 
ibrations by ISchmidt-Kalerl Jl982h . IWalbornl ll972h . and 
Garmany & Stencell dl992h ~respectively. On average, our dd yn 
lead to systematically lower My than the previous calibrations. 



icant parallax in Ivan Leeuwenl (I2007I) . has an My, s t - M VtHip p 
-1.34. 



5. Individual stars 



In the following we discuss the results for a few selected 
stars in more detail. These include the case of ( Pup, the closest 
O-star to the Sun, a case of a common origin of 4 O-stars, and 
finally a case of O-stars presently at large distance from the 
Galactic plane. 



5.1. £ Pup (HD 66811) 

( Pup is an important benchmark for the astrophysical char- 
acteristics of massive stars . Spectr o scopic ally it is classi- 



fied as 04 1 according to ISota et alj d2007l) . and its trigono 



metri c parallax from the original Hipparcos catalogue (ESA 
1997b is 2.33+0 . 51 m as. From a back-tracing of £ Pup, 



Hoogerwerf et al. (2001) found that this star had a possible en 
counter with the cluster Trumpler 10 some 2 Myr ago provided 
that its dynamical distance was dd vn = 250 ... 300 pc. This was 
inconsistent with the Hipparcos distance of 430 pc, and also 
its absolute magnitude did not agree with the cluster isochrone. 



Our b ack-tracing confirms the results from iHoogerwerf et al 
d200ll) . giving Trumpler 10 as the host and dd yn = 300 pc, 
tenc - 2.5 Myr, d enc = 0.9 pc, and pkm = 0.94. This re- 
sult is consistent with the new Hipparcos paral lax (3.00+0.1 
mas) from the re-reduction of Hipparcos data by Ivan Leeuwen 
(2007). If we literally adopt the new Hipparcos distance of 
333 pc, we get a solution with t enc = 1.8 Myr and d enc = 7.1 pc, 
which has only a slightly smaller probability (p^„ = 0.91). 

In Table Q] we summarise the distances of £ Pup from the 
various sources. Compared to the new Hipparcos benchmark, 
the old spectroscopic distance renders ( Pup 1.34 mag too 
bright in absolute magnitude. According to our dd yn it would be 
0.24 mag fainter than based on the revised Hipp arcos distance. 



From a non-LTE analysis of the spectrum, iKudritzki et al 



Table 1. £ Puppis: distances, distance moduli V - M v , and de- 
rived differences in abso lute magnitude M v with respect to the 
new Hipparcos parallax ( van Leeuwenll2007 ). 



distance, pc 


V- M v 


AM V 


sources for distance 


333 


7.60 


0.00 


van Leeuwen (2007) 


429 


8.16 


-0.56 


ESA (1997) 


615 


8.94 


-1.34 


Spectral type 04 I 


300 


7.36 


+0.24 


this paper 



= 42000 K instead of 50000 K according to spectral type 04. 
This, together with the low \ogg = 3.5, means that ( Pup is 
alr eady away from t he ZA MS . 



(1983) found that the effective temperature of f Pup is T e ff 



Kudritzki et al.l (119831) also determined the angular diame- 
ter of £ Pup to be a = 4.0 x 10~ 4 arcsec. Together with an as- 
sumed distance of (450± 200) pc, this yielded a radius of (19 + 
8)/? Q and a mass of 40 M . Using the new Hipparcos parallax, 
we find a radius of (14 ± 0A)R Q and a mass of (22.5+1.3) M . 
For its luminosity we find logL/L Q = 5.74 ± 0.02. For the error 
calculation, only the mean error of the Hipparcos parallax is 
consid ered. The new Hipparcos pa rallax rules out the scenar- 
ios by Ivan Rensbergen et al. (1996). They discussed an origin 
of £ Pup in Vela R2, which would imply a present-day dis- 
tance of the star of 700 to 800 pc. The alternative scenario 
that £ Pup originated as a field star and its runaway nature 
is due to a binary history is now also ruled out, because this 
would lead to a present-day distance between 400 and 800 pc 
(Ivan R ensberge n et al ] ll996l) . 

To summarise, £ Pup is closer, less massive, and less lumi- 
nous than previously thought. 

5.2. HD 188209, HD 189957, HD 198846, HD 201345 

Although, at present, these stars are separated by hundreds 
of parsecs, they turn out to have common "relatives". We 
show their spatial distribution in Fig. [8] where X, Y, Z 
are the Cartesian Galactic coordinates and RG(pc) is the 
Galactocentric radius of a star or a cluster. The present-day lo- 
cation is shown in the left column (f = 0), whereas the right 
column (f = -10 Myr) displays the same region 10 Myr ago. 
At that time the open c luster NGC 6871 had fo rmed (its age 
is dated logf = 6.99 in IKharchenko et al.ll2005al) . Just before, 
at te„c = 10.9 Myr, the O-star HD 188209 escaped from this 
region with a relative velocity of about 35 km/s, away from the 
Galactic plane, and towards the North pole. A few million years 
later (f = -7.7 Myr, not shown in Fig. ID, NGC 6871 passed at 
about 35 pc the region where the association Cyg OB 2 came 
into being at t « -5.5 Myr. IKharchenko et al . (2005a) give 
log t = 6.72 for Cyg OB 2. Immediately before, two O-stars 
were ejected from this region, HD 189957 (t enc =5.9 Myr) and 
HD 201345 (t enc = 5.8 Myr). HD 189957 started in the direc- 
tion to the North pole with a relative velocity of about 70 km/s 
whereas HD 201345 moved to the South pole with about the 
same relative velocity. About 3.5 Myr later (f = -2 Myr) the 
O-star HD 198846 was ejected from Cyg OB 2 with a relative 
velocity higher than 100 km/s in the direction to the South pole, 
too. 
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Fig. 8. Spatial distribution of the O-stars (1: HD 188209, 2: 
HD 189957, 3: HD 201345, 4: HD 198846), of the open cluster 
NGC 6871, and of the association Cyg OB 2. The upper panels 
are for the (X, F)-plane, the lower panels show the distribution 
in the (RG,Z)-plane, where RG is the Galactocentric radius. 
Asterisks mark the location of stars, the circles represent the 
clusters; th eir radii are defined as t he present-day tidal radii de- 
termined in lPiskunov et alJ (I2008al) . The left columns show the 
location at present; the second, third, and forth columns show 
the locations 2 Myr, 5.5 Myr, and 10 Myr ago, respectively. 
The dashed curves in these panels delineate the stellar orbits 
computed with the back-tracing procedure from present to the 
corresponding time t indicated on the top. Note that Cyg OB 2 
is younger than 10 Myr, and therefore does not appear in the 
right-most column. 



As we note in Sec. 3 and show in Fig. [2] there is, al- 
though with a slightly smaller kinematic probability, another 
acceptable solution for HD 201345 suggesting an ejection from 
NGC 6871 about 6.5 Myr ago. Also for HD 189957, a second 
acceptable solution turned out to be possible (an ejection from 
NGC 6871 with the parameters d dyn = 1964 pc, t enc =7.1 Myr, 
dene - 3.5 pc, pkin = 0.84). Since the time of their closest 
approach with NGC 6871 is comparable with the time when 
NGC 687 1 passed the birth place of Cyg OB 2, one may suspect 
that this event was the reason for the ejection of HD 201345 
and HD 189957 from a region in between NGC 6871 and 
Cyg OB 2. 

5.3. O-stars at large distances from the Galactic plane 

The majority of O-stars in our sample is located within 200 pc 
from the Galactic plane. However, there are five stars with 
a present-day location at |Z| > 400 pc. These stars either 
had relative velocities at the moment of encounter (ejection) 
larger than 100 km/s (HD 116852, HD 157857) or they left 
their cluster counterparts more than 9 Myr ago (HD 14633, 
HD 105056, HD 175754). Below, we consider two examples, 
HD 116852 and HD 14633, which have the largest distances 



from the Galactic plane. These cases are also interesting since 
their counterparts were clusters hosting other field O-stars, too. 

5.3.1. HD 116852 and Trumpler 14 

Judging from the spectroscopic distance and high galactic 
latitude of HD 116852, one can expect that this star is lo- 
cated at a relatively large distance from the Galactic plane. 
Indeed, from back-tracing we obtain Z =s -690 pc for the 
present-day location of HD 116852. This star left its birth- 
place in a protocluster of Trumpler 14 with a relative veloc- 
ity of » 180 k m/s about 6 Myr ago, i.e. j ust before the clus- 
ter was formed. Kharchenko et al. ( 2005ab list log? = 6.67 for 
Trumpler 14. The location of the star and of the protocluster 
at that time was at Z « -50 pc. During the next 6 Myr, the 
cluster moved by only 60 pc, whereas HD 116852 did more 
than 1 kpc. Merely 1.5 Myr after its birth, Trumpler 14 lost an- 
other O-star, HD 93652, and about 1.5 Myr later, HD 91651 
and HD 305539 left the cluster, too. Presently, HD 93652, 
HD 91651, and HD 305539 are located about 40 pc, 75 pc, 
and 30 pc away from Trumpler 14, respectively. According to 
Mel'nik & Efremovl(ll995h . there is an OB-association Car 1 F 
at about 65 pc from Trumpler 14. Probably, this neighbourhood 
has had impact on the fate of Trumpler 14. 

5.3.2. HD 14633 and ASCC 8 

Though presently HD 14633 is at Z * -670 pc, about 9 Myr 
ago this star was located near the Galactic plane between the 
clusters NGC 869 (h Per) and NGC 884 Of Per), i.e. within a 
region wellknown as the then very active star formation region 
Per OB 1. About 3.5 Myr later this place was also the birth- 
place of an open cluster, ASCC 8 (logf = 6.76). In Fig. 9 we 
show the location of these three clusters (NGC 869, NGC 884, 
ASCC 8) as it is now and as it was in the past. 



Acc ording to the tidal radii determined in Pis kunov et al 



(2008a), NGC 869 and NGC 884 are the largest clusters in the 
whole COCD cluster sample. Moreover, their tidal spheres are 
overlapping, and some concentration of stars above the back- 
ground is observed in this overlapping zone. Based on com- 
mon proper motions, an open c luster (ASCC 8) was identified 
here dKharchen ko et al. 2005b), with an age of about 5.8 Myr. 
Obviously, such a location of a cluster does not present a good 
chance for a long life. From the back-tracing computations, 
three O-stars marked as field in the GOSV2 came from the 
region occupied now by ASCC 8: HD 14947 - 1.2 Myr ago, 
HD 17603 - 7.3 Myr ago, and HD 14633 - 9.2 Myr ago. Two 
other field stars, HD 12993 and HD 13022, had their origins in 
NGC 884 and NGC 869, respectively. Note that the final solu- 
tion for a given star and its host is selected from all acceptable 
solutions (de„c < 10 pc) to have the smallest encounter distance 
dene at the highest kinematic probability pu„. This approach is 
justified if possible hosts (i.e. open clusters) are separated by 
distances which are considerably larger than their sizes. This 
is the usual case. However, the example considered presents 
an exception from the rule. Here the distances between clus- 
ter centres are comparable to their tidal radii. Therefore, so- 
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Fig. 9. Spatial distribution of the field O-stars (1: HD 14633, 
2: HD 17603, 3: HD 14947, 4: HD 12993, 5: HD 13022), and 
of the open clusters NGC 869, NGC 884, ASCC 8. The left 
panel is for the {X, 5 / )-plane, whereas right panel shows the dis- 
tribution in the {RG, Z)-plane, RG is the Galactocentric radius. 
Asterisks mark the present-day location of the stars, the circles 
are for the clusters: solid curves for the present-day location 
and dash-dot-dot-curves for the past. Their ra dii are defined 
by the present-day tidal radii determined in Piskunov et al.1 
(2008a). The dashed and dotted curves show the orbits of the 
stars and clusters, respectively, computed with the back-tracing 
procedure. The crosses mark the places where the stars encoun- 
tered their hosts. Negative numbers in the left panel express the 
time in Myr when the corresponding events happened (in order 
not to overload the figure, the numbers are omitted in the right 
panel). 

lutions with d enc > 10 pc can be considered, too. For exam- 
ple, HD 12993 could have been ejected from NGC 869 about 
0.8 Myr ago when the star was « 50 pc away from the cluster 
centre. In this special case, we prefer a more general statement: 
the observational data for these five field O-stars are consistent 
with the claim that they originated in Per OB 1 . 

6. Stars without solutions 

For 73 out of 105 O-stars with assignments as field stars 
in the GOSV2 we could trace back their origins and found the 
corresponding host clusters orprotoclusters. Twelve stars could 
not be treated by our method, because no radial velocities were 
available for them. For another twenty, we could not find a so- 
lution, i.e. they could not be associated with any of the clusters 
of our sample in the past 10 million years. This negative find- 
ing, however, should not be interpreted as a proof that these 
O-stars have formed in isolation. 

We could not find a parameter or parameters which dis- 
tinguish these stars from the other field O-stars which had an 
acceptable solution. Both groups cover a similar range of ap- 
parent magnitudes, and they are comparable in the distribution 
of the mean errors of their kinematic components. Of course, 
one cannot exclude that a few of them have a true distance 
modulus differing from the spectroscopic estimates by more 
than two magnitudes, and/or their true velocity components 
are outside the intervals checked with the backward procedure. 
Nevertheless, a more important aspect seems to be that our 
sample of potential host clusters is not complete for associa- 



tions and very young (embedded) clusters, as we stressed in 
Sec. 2.1. Provided that this assumption is correct, one would 
expect about 20% of field O-stars having their origins in these 
kinds of objects. To answer the question with certainty, how- 
ever, more accurate data on the distances and kinematics of 
associations and embedded clusters are required. 



7. Summary 



In this paper we have followed the dynamical history of 
O-stars that left the groups where they originated in. We leave 
open the physical mechanisms which are behind these events. 

For 73 out of 93 O-stars considered, we found accept- 
able solutions indicating that the present-day data are consis- 
tent with the assumption that these O-stars were ejected from 
young open clusters or protoclusters during the past 10 Myr. 
The GOSV2 catalogue counts 370 O-stars, for 265 of which 
the origin is given in that catalogue. We were able to add 73 
more cases to the list. For 32 stars (or 9%) we could not prove 
the origin in groups. 

In this paper we have dealt with aspects of the early phases 
in the life of open clusters. O-stars are best suited as tracers 
of this early-phase evolution because of their short life-time. In 
our sample of 16 1 young open clusters (age < 30 Myr) from the 
COCD there are 55 (or » 35%) that have O-stars as members 
( Kharchenko et al. 2005 allbh . 23 of these have already lost one 
or more O-stars in their history. Another 24 (or 15%) of the 
COCD clusters had relations to O-stars in the past 10 Myrs, but 
do not contain O-stars at present. For 82 (5 1 %) young clusters 
we cannot prove a relation to presently living O-stars. Either 
their most massive member is a main sequence star of spectral 
type later than O, or it is a former O-star which has already 
evolved. Of the 47 clusters that have lost at least one O-star, we 
find 14 that are so young that O-star and cluster should already 
have separated in the protocluster phase. 

Summing up the statistics above, the following picture 
emerges. Fifty percent of the clusters being able to survive 
the infant-mortality phase are so massive that they contain or 
contained O-stars. These O-stars have not been able to de- 
stroy the clus t er. Th is, in parts, answers the question asked by 
Lada & Lada d2003h : Do the progenitors of bound open clus- 
ters ever contain O-stars? Yes, they did. 

The fact, that we could not trace back 9% of all O-stars 
from the GOSV2, does not necessarily mean that "isolated 
O-star formation" is possible. Their known astrophysical data 
(distance, velocities) may be incorrect or our list of possible 
host candidates may be incomplete. On the other hand, we can 
interpret our result - no solution for 20 out of 93 stars - as fol- 
lows: there is an upper bound of slightly more than 20% of O- 
stars which could have destroyed their family of brothers and 
sisters with which they may have fo rmed to gethe r. 

It has been shown by lPiskunov et a l. (2008b) that classical 
(gravitationally bound) open clusters in the Milky Way evolve 
due to stellar and dynamical evolution as well as due to external 
perturbations. They are losing stellar mass during their live- 
time at an average rate of 3 to 14 M /Myr. In this paper we 
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determined the mass loss rate of young open clusters due to O- 
stars alone to be 3 to 5 M /Myr in the first few million years of 
their existence. 

As a by-product, we find new distances and absolute mag- 
nitudes for 73 O-stars. These indicate that the calibration of 
absolute magnitudes of O-stars should be revised. Their ab- 
solute magnitudes are systematically fain ter by about 0.3 to 
0.8 m a g compared to t he calibrations by Garmany & Stencell 
dl992hjwalbornl dl972l) . or lSchmidt-Kalerl dl982b . This would 
be consistent with the re-calibration of the a bsolute magnitudes 
of early B-type stars by iKaltcheva (1999) using Hipparcos 
trigonometric parallaxes. 

We have also shown that ( Pup, the closest O-star from the 
Sun, left the young open cluster Trumpler 10 some 2.5 Myrs 
ago. Its present-day distance from the Su n of 300 pc is compati - 
ble with the new Hipparcos distance from lvan Leeuwen1d2007l) . 
This implies a radius of 14 R Q , a mass of 22.5 M Q , and a lumi- 
nosity of logL/L of 5.74 for £ Pup, i.e. the values are consid- 
erably smaller than assumed before. 

Acknowledgements. We are grateful to Nina Kharchenko, Anatoly 
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formation and young open clusters. 
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HD 


Sp. type 


Cluster name 


log t oc l 


ddyn 


tenc 


dene 


AVel enc 


My 


min, 


max] 


Pkin 


comments 










[pc] 


[Myr] 


[pc] 


[km/s] 




[mag 


I 






1337 


O9.0 III 


NGC 957 


6.84 


1911 


8.5 


3.9 


75.7 


-6.0 


[-6.4, 


-5.8] 


0.54 




10125 


09.7 II 


IC 1590 


6.84 


2300 


4.9 


6.1 


78.5 


-5.3 


[-5.6, 


-5.0] 


0.81 




12323 


O9.0 V 


ASCC9 


6.79 


2890 


6.7 


4.6 


46.8 


-4.0 


[-4.1, 


-3.5] 


0.97 




12993 


06.5 V 


NGC 884 


7.10 


2074 


4.3 


0.3 


61.2 


-4.0 


[-4.2, 


-3.6] 


0.71 




13022 


09.5 II-III 


NGC 869 


7.28 


1948 


8.3 


2.6 


17.2 


-4.3 


[-4.7, 


-4.1] 


0.88 




13745 


09.7 II 


NGC 663 


7.14 


2035 


4.7 


2.5 


50.5 


-4.9 


[-5.2, 


-4.6] 


0.80 




14434 


05.5 V 


ASCC9 


6.79 


3062 


5.9 


5.5 


34.0 


-5.2 


[-5.5, 


-4.9] 


0.81 




14442 


O5.0 V 


ASCC9 


6.79 


2912 


3.9 


3.4 


40.9 


-5.2 


[-5.4, 


-4.8] 


0.94 




14633 


O8.0 V 


ASCC 8 


6.76 


2151 


9.2 


3.8 


76.5 


-4.5 


[-4.6, 


-4.0] 


0.89 


(1) 


14947 


O5.0I 


ASCC8 


6.76 


2177 


1.2 


8.8 


32.1 


-5.8 


[-6.1, 


-5.5] 


0.77 




15137 


09.5 II-III 


NGC 957 


6.84 


1784 


5.8 


2.4 


27.9 


-4.3 


[-4.5, 


-3.9] 


0.92 




16691 


O4.0I 


ASCC 9 


6.79 


2901 


1.1 


6.8 


43.3 


-5.8 


[-6.1, 


-5.5] 


0.55 




17603 


07.5 lb 


ASCC 8 


6.76 


2272 


7.3 


1.2 


26.5 


-6.1 


[-6.3, 


-5.7] 


0.97 




39680 


O6.0 V 


NGC 2169 


6.89 


1354 


9.1 


5.5 


35.2 


-3.6 


[-3.8, 


-3.2] 


0.84 




41997 


07.5 V 


Collinder 89 


7.50 


723 


2.8 


2.4 


44.4 


-2.7 


[-3.0, 


-2.5] 


0.63 




44811 


O7.0 V 


NGC 2129 


7.08 


1582 


9.5 


8.9 


18.3 


-3.7 


[-4.0, 


-3.4] 


0.51 




52266 


O9.0 IV 


Collinder 106 


6.74 


1388 


10.4 


1.9 


34.3 


-4.3 


[-4.4, 


-3.8] 


0.55 




52533 


09.5 V 


Collinder 106 


6.74 


1831 


10.0 


3.8 


34.9 


-4.2 


[-4.4, 


-3.7] 


0.86 




57236 


O8.0 V 


NGC 2414 


6.94 


2432 


9.1 


5.4 


62.3 


-4.6 


[-4.8, 


-4.2] 


0.84 




66811 


O4.0I 


Trumpler 10 


6.95 


297 


2.5 


1.4 


51.0 


-5.2 


[-5.3, 


-4.7] 


0.94 


(1) 


69464 


06.5 lb 


ASCC 45 


7.12 


3001 


0.1 


0.8 


11.4 


-5.2 


[-5.6, 


-4.9] 


0.78 


(2) 


75222 


09.7 lab 


Collinder 205 


7.03 


1735 


6.6 


1.1 


62.8 


-5.6 


[-6.0, 


-5.4] 


0.97 


(1) 


76968 


09.7 lb 


ASCC 45 


7.12 


2682 


10.6 


4.1 


81.2 


-5.9 


[-6.5, 


-5.8] 


0.86 




89137 


09.5 III 


Loden 306 


6.76 


2026 


9.1 


3.1 


45.6 


-4.3 


[-4.4, 


-3.8] 


0.91 




91651 


O9.0 V 


Trumpler 14 


6.67 


2720 


1.6 


2.3 


46.3 


-4.2 


[-4.3, 


-3.7] 


0.91 




93632 


O5.0 III 


Trumpler 14 


6.67 


2731 


3.1 


6.0 


12.3 


-5.6 


[-5.8, 


-5.1] 


0.74 




94963 


06.5 III 


IC 2581 


7.22 


2504 


8.7 


6.5 


27.8 


-5.1 


[-5.6, 


-5.0] 


0.81 




96917 


08.5 lb 


Collinder 228 


6.68 


1991 


3.1 


1.9 


55.5 


-5.5 


[-5.7, 


-5.1] 


0.84 




96946 


O6.0 V 


vdBergh-Hagen 121 


6.64 


2738 


4.1 


7.8 


49.2 


-5.3 


[-5.4, 


-4.8] 


0.70 




97848 


O8.0 V 


NGC 3324 


6.72 


2329 


4.0 


8.8 


48.0 


-4.0 


[-4.1, 


-3.5] 


0.90 




104565 


09.7 la 


ASCC 75 


6.65 


2839 


4.7 


3.2 


148.0 


-5.1 


[-5.2, 


-4.5] 


0.72 




104649 


09.5 V 


NGC 3572 


6.88 


1948 


5.8 


2.2 


40.2 


-4.4 


[-4.6, 


-4.0] 


0.93 




105056 


09.7 la 


Ruprecht 94 


7.19 


3304 


9.0 


3.8 


47.4 


-5.8 


[-6.4, 


-5.7] 


0.64 


(1) 


105627 


O9.0 II-III 


Loden 821 


7.29 


2885 


6.5 


7.5 


70.8 


-4.8 


[-5.3, 


-4.7] 


0.68 




1 12244 


08.5 lab 


Feinstein 1 


6.97 


1391 


8.8 


3.2 


45.3 


-6.0 


[-6.4, 


-5.8] 


0.83 




116852 


O9.0 III 


Trumpler 14 


6.67 


2475 


6.0 


2.7 


179.6 


-4.1 


[-4.2, 


-3.6] 


0.67 


(1) 


117856 


09.5 III 


Loden 694 


7.38 


1750 


10.0 


5.2 


17.4 


-5.0 


[-5.5, 


-4.9] 


0.51 




120521 


O8.0 lb 


Loden 821 


7.29 


2890 


7.8 


1.8 


26.5 


-5.0 


[-5.5, 


-4.8] 


0.80 




120678 


O8.0 III 


Loden 821 


7.29 


2864 


5.9 


1.5 


39.8 


-5.2 


[-5.8, 


-5.1] 


0.84 




123008 


09.7 lb 


Loden 821 


7.29 


2875 


6.9 


1.4 


51.0 


-5.1 


[-5.6, 


-5.0] 


0.66 




123056 


09.5 V 


Loden 694 


7.38 


1577 


10.1 


3.1 


28.5 


-3.9 


[-4.3, 


-3.7] 


0.83 




125206 


09.5 IV 


NGC 5606 


6.84 


2054 


7.2 


3.6 


41.1 


-5.2 


[-5.5, 


-4.9] 


0.68 




135240 


07.5 III 


ASCC 79 


6.86 


809 


2.2 


0.6 


5.2 


-5.0 


[-5.3, 


-4.7] 


0.73 


(3) 


135591 


07.5 III 


ASCC 79 


6.86 


796 


0.3 


1.4 


17.1 


-4.6 


[-4.8, 


-4.2] 


0.74 


(3) 


148546 


O9.0 la 


ASCC 88 


7.17 


1453 


8.8 


5.7 


55.4 


-4.7 


[-5.0, 


-4.4] 


0.73 


(1) 


153426 


O9.0 II-III 


Hogg 22 


6.70 


1694 


8.1 


9.9 


55.4 


-4.9 


[-5.1, 


-4.5] 


0.91 




153919 


06.5 la 


NGC 6231 


6.81 


1034 


5.1 


1.1 


48.8 


-5.1 


[-5.4, 


-4.7] 


0.82 




154368 


09.5 lab 


Sco OB4 


6.82 


1082 


4.7 


5.3 


11.2 


-6.3 


[-6.6, 


-6.0] 


0.65 




154643 


09.5 III 


ASCC 88 


7.17 


1908 


1.4 


2.3 


20.3 


-5.6 


[-6.1, 


-5.5] 


0.85 




154811 


09.7 lab 


vdBergh-Hagen 205 


7.12 


2019 


8.2 


2.7 


35.1 


-6.2 


[-6.8, 


-6.2] 


0.96 




156212 


09.7 lab 


Trumpler 28 


6.89 


992 


7.0 


4.0 


53.5 


-4.5 


[-4.7, 


-4.1] 


0.93 


(1) 


157857 


06.5 III 


NGC 6611 


6.72 


1854 


3.8 


2.6 


113.5 


-4.9 


[-5.1, 


-4.5] 


0.68 


(1) 


158186 


09.5 V 


Sco OB4 


6.82 


1254 


5.6 


1.2 


30.3 


-4.4 


[-4.6, 


-4.0] 


0.82 




161853 


O8.0 V 


Trumpler 28 


6.89 


1241 


5.0 


2.6 


25.4 


-4.1 


[-4.3, 


-3.7] 


0.68 
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166734 


O7.0 lb 


NGC 6604 


6.64 


1661 


3.0 


2.5 


22.4 


-6.7 


;-7.o, 


-6.4] 


0.65 




169515 


09.7 lb 


Markarian 38 


6.95 


1794 


8.9 


4.3 


37.7 


-6.1 


>6.5, 


-5.9] 


0.54 




169582 


O6.0I 


NGC 6604 


6.64 


1565 


5.5 


3.2 


28.5 


-4.9 


>5.0, 


-4.4] 


0.57 




171589 


O7.0 II 


NGC 6618 


6.78 


1934 


3.3 


1.5 


63.2 


-4.9 


;-5.i, 


-4.5] 


0.96 




175754 


O8.0 II 


ASCC 93 


7.22 


2586 


10.3 


4.9 


64.8 


-5.3 


>5.8, 


-5.2] 


0.47 




175876 


06.5 III 


NGC 6618 


6.78 


1922 


5.5 


4.8 


61.7 


-5.0 


>5.2, 


-4.5] 


0.88 




188001 


07.5 la 


Roslund 2 


6.89 


2022 


2.7 


2.5 


70.3 


-6.0 


;-6.4, 


-5.7] 


0.79 


(1) 


188209 


09.5 lab 


NGC 6871 


6.99 


1705 


10.9 


7.3 


33.7 


-5.9 


>6.4, 


-5.7] 


0.57 




189957 


09.5 III 


Cyg OB2 


6.72 


1858 


5.9 


3.1 


67.7 


-4.4 


>4.6, 


-4.0] 


0.87 


(1) 


192281 


05. V 


NGC 6913 


7.12 


1081 


2.4 


1.2 


35.5 


-4.5 


-4.8, 


-4 21 


0.80 




195592 


09.7 la 


NGC 6913 


7.12 


1132 


2.9 


6.3 


40.3 


-6.3 


;-6.9, 


-6.3] 


0.79 




198846 


O9.0 V 


Cyg OB2 


6.72 


1427 


1.9 


3.7 


106.6 


-4.1 


;-4.3, 


-3.7] 


0.81 


(1) 


201345 


O9.0 V 


Cyg OB2 


6.72 


1739 


5.8 


1.5 


65.6 


-4.0 


;-4.i, 


-3.5] 


0.99 




203064 


07.5 III 


Collinder419 


6.85 


631 


6.5 


1.3 


26.9 


-4.8 


;-5.o, 


-4.4] 


0.87 


(1) 


218915 


09.5 lab 


ASCC 120 


7.08 


2369 


5.8 


4.3 


54.7 


-5.4 


;-5.8, 


-5.1] 


0.72 


(1) 


303492 


O9.0 la 


Loden 153 


6.74 


2714 


1.5 


2.6 


79.3 


-5.8 


;-6.i, 


-5.5] 


0.71 




305523 


O9.0 II 


ASCC 65 


7.09 


3532 


7.5 


6.3 


27.8 


-5.3 


;-5.7, 


-5.1] 


0.81 




305532 


O6.0 V 


ASCC 65 


7.09 


3448 


3.6 


6.7 


58.0 


-4.3 


;-4.6, 


-4.0] 


0.80 




305539 


O7.0V 


Trumpler 14 


6.67 


2752 


1.4 


1.3 


21.3 


-4.1 


;-4.2, 


-3.6] 


0.83 





Comments: 



(1) - runaway O-stars according to lSota et al. I d2007|) 



(2) - member of the ASCC 45 cluster according to Kharchenko et al 



(3) - member of the ASCC 79 cluster according to Kharch enko et al 



|Khi 
iKhi 



(2005b) 
(2005b) 



